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ABSTRACT 

Electroproduction of pions on the nucleon near the threshold is analyzed within the frame- 
work of baryon chiral perturbation theory. We give a thorough discussion of the low-energy 
theorems related to charged and neutral electropionproduction. It is shown how the axial 
radius of the nucleon can be related to the S-wave multipoles E^j^ and L^j^ . The chi- 
ral perturbation theory calculations of the 7*p — *• n^p reaction are found to be in good 
agreement with the recent near threshold data. We also discuss the influence of some 
isospin-breaking effects in this channel. For future experimental tests of the underlying 
chiral dynamics, extensive predictions of differential cross sections and multipole ampli- 
tudes are presented. 
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I. INTRODUCTION AND SUMMARY 

Threshold pion photo- and electroproduction has received renewed interest over the 
past few years from both the experimental and the theoretical side. Extensive data were 
obtained for the process 7 + p — 7r° + p. These data spurred numerous theoretical investi- 
gations concerning the low-energy theorem for the electric dipole amplitude £^0+ and also 
a critical reexamination of the data was performed. In the framework of QCD these topics 
were addressed in big detail in Ref.[l] (which from now on will be referred to as I). Using 
virtual photons as probes one can get more detailed information about the structure of the 
nucleon in the non-perturbative regime of QCD due to the longitudinal coupling of the 
virtual photon to the nucleon spin. For example, charged pion electroproduction encodes 
information about fundamental quantities like e.g. the axial form factor of the nucleon. 
Many of these topics are addressed in the monograph by De Alfaro et al. [2]. Also, new 
experimental information has recently become available for the process 7* + p — > tt^ + p, 
where 7* denotes the virtual photon [3] . This experiment is a major step beyond previous 
measurements which were characterized by poor energy resolution and did not come close 
enough to the production threshold. Therefore, the results of the old measurements were 
dominated by the P-wave Mi+ multipole. In contrast, the data presented in [3] were ob- 
tained in the energy range of to 2.5 MeV above threshold and at photon four-momenta 
squared of k'^ = -0.04 to -0.1 GeV^ (A;^ < in the physical region of the electroproduc- 
tion reaction). More accurate data also for the production of charged pions are expected 
to come from NIKHEF, MAMI at Mainz and MIT Bates in the near future. 

In this paper, we will be concerned with a systematic analysis of the processes 7*+p — > 
TT'^ + n, 7* + n — > 7r~ + p and 7* + p — > 7r° + p in the threshold region making use of 
baryon chiral perturbation theory (CHPT). In general, CHPT allows one to systematically 
investigate the strictures of the spontaneously broken chiral symmetry of QCD. It is based 
on the observation that in the three flavor sector of QCD, the quark masses are small and 
that the theory in the limit of vanishing quark masses admits an exact chiral symmetry. 
The latter is dynamically broken which leads to the appearance of massless pseudoscalar 
excitations, the Goldstone bosons. In the real world, the quark masses are not exactly zero 
and thus the Goldstone bosons acquire a small mass. The interaction of these particles with 
each other and matter fields like e.g. the nucleons are weak at low energies as mandated by 
Goldstone's theorem. This fact is at the heart of CHPT which amounts to a systematic and 
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simultaneous expansion of the QCD Green functions in small momenta and quark masses. 
To perturbatively restore unitarity it is mandatory to consider pion loop diagrams. In what 
follows, we will work in the one-loop approximation which has been shown to be of sufficient 
accuracy for many threshold phenomena. For a general review see Ref.[4]. At next-to- 
leading order, one has not only pion loop graphs but also local contact terms. The latter are 
accompanied by a priori unknown coefficients, the so-called low energy constants. These 
can either be determined from phenomenology [5] or estimated from resonance exchange 
[6]. For the case under consideration, A (1232) and vector meson exchanges, the knowledge 
of some nucleon electromagnetic radii as well as the pion charge radius allow us to pin 
down all low energy constants appearing here. 

The starting point of our investigation will be flavour SU(2) with equal masses and we 
work to first order in the electromagnetic coupling constant. As a consequence the one- 
loop approximation does not include isospin-breaking effects as they are revealed in the 
difference of the proton and neutron or the charged and neutral pion masses. To minimally 
account for these effects we will also present calculations in which the respective particles 
have their physical masses. For the neutral pion photoproduction it was argued in [7] that 
indeed the pion mass difference is the most important effect of the isospin breaking. In 
that reference it was also shown that this procedure leads to a much improved description 
of the data for Eqj^^ in the threshold region. We should stress, however, that ultimately 
a complete calculation including effects of higher loops and of higher order in the fine 
structure constant should be performed. This is beyond the scope of the present paper. 

The pertinent results of our investigation can be summarized as follows: 

(i) To one-loop order and to first order in the electromagnetic coupling, there are 66 
topologically different Feynman diagrams including pion loops. These diagrams can 
be divided into three separately gauge invariant subsets (the first class was already 
discussed in ref.[14]). After mass and coupling constant renormalization, the corre- 
sponding invariant amplitudes A\'^^ {i = 1, 2, 3, 4, 5, 6; a = +,0,—) contain further 
UV divergences which can be absorbed by two low-energy constants related to the 
isovector charge radius of the nucleon and the electromagnetic radius of the pion. In 
addition, there are six finite contact terms of order q^. Four of the unknown coeffi- 
cients were already determined in our study of threshold pion photoproduction [1] and 
the others can be fixed from the axial and the isoscalar charge radii of the nucleon. 
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(ii) We have discussed in some detail the low-energy theorems (LETs) for the electric 
dipole amplitude Eq^ and the longitudinal amplitude Lq^ (a = +, 0, — ). For neutral 
pion electroproduction, virtual pions generate infrared singularities which modify the 
familiar form of the energy expansion [2,17]. In QCD, the expansion for Eq^ and Lq^ 
in powers of = M-^/m (the ratio of the pion to the nucleon mass) and v = k^/m^ 
(with the four-momentum of the virtual photon) is given in eq.(5.1). For 7*p 7r°p 
and 7*n — > 7r°n, these agree with the result given in ref.[14]. For the (— ) amplitudes 
which are probed in charged pion electroproduction, we have found a modification of 
the LET due to Nambu et al. [18]. This was already discussed in ref.[19] where it was 
shown that an additional pion loop effect modifies the relation between dEQ^{k'^) /dk'^ 
at k^ = and the nucleon axial radius < r\ >^/^ (already in the chiral limit). In fact, 
this correction leads to a better agreement between the axial radius determined from 
(anti) neutrino scattering [13] and from pion electroproduction [31]. For the planned 
high precision charged electroproduction experiments it is mandatory to include this 
effect in the analysis of the data. We have also worked out the terms of order /j,^ and 
fi^lnfi^ from the one-loop graphs for Lq_^{ii,0) and shown that the expansion in 
converges slowly. This agrees with our findings for photoproduction reported in ref.[l]. 

(iii) We have confronted the chiral prediction with the recent very accurate NIKHEF data 

for 7*p n^p [3]. As already shown in ref.[30], one-loop effects are necessary to 
understand the /c^-dependence of the S-wave cross section for photon four-momenta 
|/c^| < 0.1 GeV^. The theoretical prediction for Lo+ at the photon point, |ivo+P = 
0.2 //6 is in fair agreement with the result of ref.[3], |-^o+P = 0.13 ± 0.05 ^b. 

(iv) The calculated multipole amplitudes are presented. For future experimental tests of 
our predictions, extensive results are also given for the transverse (T), longitudinal 
(L), interference (/) and polarization (P) differential cross sections daT,L,i,p/dO, and 
are shown for all three channels, i.e. 7*p tt'^ti, j*n — > tt~p and 7*p — > ir'^p. To 
see the dynamical content of our CHPT calculations, these results are compared with 
the traditional calculations using pseudo-vector pion coupling with the nucleon. Our 
results strongly suggest that to test the dynamical content of the chiral perturbation 
theory prediction, it is mandatory to perform a transverse/longitudinal separation of 
the data. We have also shown that for photon four-momenta \k'^\ > 0.05 GeV^, the 
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loop corrections become large and thus it would be preferable to have experiments at 
lower photon four-momenta. 

(v) For neutral electropionproduction of the proton, we have investigated the dominant 
effect of isospin breaking, namely the charged to neutral pion mass difference in the 
various loop functions. We have given a gauge invariant prescription to implement 
this effect. The resulting S-wave cross section is marginally different from the one 
presented in ref.[30]. Only at k'^ = 0, isospin breaking seems to be of importance (due 
to the dominance of the class I diagrams). We get an improved prediction for Lo+ 
at the photon point, |i^o+P = 0.f7/u6. At finite fc^, we could only find a substantial 
effect on dap/dVl (at /c^ = —0.04 GeV^). The largest uncertainty therefore resides in 
the knowledge of the low-energy constant di (see ref.[l]). 

The paper is organized as follows. In section 2 we discuss some formal aspects of 
pion electroproduction to fix our notation. The effective Lagrangian which will be used is 
given in section 3. We only exhibit the new terms related to the electroproduction reaction 
and refer the reader to I for more detailed discussions. Section 4 contains the invariant 
amplitudes for the reaction 7*p TT^p. Pion electroproduction low energy theorems are 
derived and discussed in section 5. We also review critically previous attempts to formulate 
these low energy theorems. The numerical results are presented in section 6 together 
with some proposals for future experiments. Various technicalities are relegated to the 
appendices. 

II. PION ELECTROPRODUCTION: SOME FORMAL ASPECTS 

In this section, we define the current matrix elements for the ^ ttN process 

in terms of invariant amplitudes which will then be subject to the chiral expansion in 
section IV. The formulas used in our calculations of the N{e, e'7r)N cross sections from 
the current matrix elements can be found in Refs.[8,9] and are therefore omitted here. To 
see the dynamics at threshold, the multipole decomposition of the current matrix elements 
at threshold is presented explicitly in this section. The formulae needed for calculating 
the current matrix elements from the multipole amplitudes tabulated in tables la,b,c are 
given in Appendix C. This will allow the readers to calculate the CHPT predictions of 
(e, e'n) spin observables accessible to possible future experiments by using the well known 
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formulae such as given in Ref.[9]. The reader famihar with these topics is invited to skip 
this section. 

II.l. BASIC CONSIDERATIONS 

Consider the process 7*(/c) + N{pi) 7r"(g) + N{p2), with N denoting a nucleon 
(proton or neutron), tt" a pion with an isospin index a and 7* the virtual photon with 
k'^ < 0. The case of photoproduction, k'^ = 0, was discussed in detail in I. The pertinent 
Mandelstam variables are s = {pi + k)'^, t = {pi — p^)^ and u = {pi — q)'^ subject to the 
constraint s + t + u = 2m^ + + k^. Here, m and denote the nucleon and pion 
mass, in order. To first order in the electromagnetic coupling, e, to which we will work, 
the transition current matrix element is given as 

J^{p2,S2;q,a\pi,si;k) = 

(2.1) 

iu2l5{l'"¥Bi + 2P^52 + 2qi'Bs + 2k'' + -fi'B^ + P'^^Be + k^'lflB-^ + q''VBs]u^ 

with Sj the spin index of nucleon j and P = {pi +p2)/2. The amplitudes Bi{s, u) have the 
conventional isospin decomposition (in the isospin symmetric limit), 

Bi{s,u) = B^^\s,u) 5az + B\-\s,u) l[r„, T3] + ^(s, w) . (2.2) 

Not all of the eight -Bi(s, u) are independent. Indeed, gauge invariance k^J^ = requires 
2k'^ [BUs, u) + 2^4(5, u)] + {s-u) B2(s, u) + 2(s + u- 2rr?) B^is, m) = 

(2.3) 

AB^{s, u) + {s- u) Bq{s, u) + Ak^ B'j{s, u) + 2{s + u- 2m^) Bs{s, tt) = 

so that one can express the transition current matrix element in terms of six independent 
invariant functions, conventionally denoted by ylj(s,tt), (z = 1, 6), 

6 

Ji" = iu2i^Y,^tMs:u)m (2.4) 

with 









M^2 
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q^k\ 
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The Ai{s,u) are related to the Bi{s,u) via 

2B2 



Ai = Bi- mBe, A 



2 



As = -Bs, A^ = ~Bg, A^^Br, (2.6) 



A.= 



. + «-2mn ' ^ ' ^ 2(t - M2) / = " ^^^/ • 



For A5 we have given two equivalent forms to make clear that the possible zeros of s + 
u — 2m^ in the physical region do not lead to a pole (infinity) of the current transition 
matrix element in the physical region. In what follows, wc will use both sets of invariant 
functions (sec section IV). Under {s ^ u) crossing the amplitudes A^^2'^^ , A3 q, b[\'^q , 
b[-^ + 2Bi-\ Bi;,],^, are even, while 4+^% A^^, sJ+'°) + 2B^+''\ Bi%%, Bi^,], are 
odd. 

In terms of the isospin components, the physical channels under consideration are 

J,(7*p- A) = JW + J(+) 
M^n ^ 7r«n) = - 
Ml*p^n+n) = V2[jl°^ + 4-^] 
J^(7*n^7r-p) = V2[j(«)-j(-)]. 

Having constructed the current transition matrix element it is then straightforward to 
calculate observables. The pertinent kinematics and definitions are outlined in refs.[8,9] to 
which the reader is referred for details (see also Berends et al. [10]). 

II.2. MULTIPOLE DECOMPOSITION AT THRESHOLD 

For the discussion of the low energy theorems in section V, we have to spell out 
the corresponding multipole decomposition of the transition current matrix element at 
threshold. In the 7*A^ center of mass system at threshold i.e. = (Mtt, 0,0,0) one can 
express the current matrix element in terms of the two 5'-wave multipole amplitudes, called 
£^0+ and Lq+, 



J = 47ri(l + //) X/{-Eo+(/i, [Lo+(//, i^) - -Eo+(//, u)]ka- k}xi 
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(2.8) 



with Xij two component Pauli-spinors for the nucleon. For the later discussion we have 
introduced the dimensionless quantities 

11= , J^=—- (2-9) 

m 

The multipole £^0+ characterizes the transverse and Lo+ the longitudinal coupling of the 
virtual photon to the nucleon spin. Alternatively to Lq+, we will also use the scalar 
multipole ^0+ defined via: 

— * 

So+is,e)=^-^Lo+{s,k'') (2.10) 

Kq 

At threshold, we can express £^0+ and Lo+ through the invariant amplitudes Ai{s,u) 
via (suppressing the isospin indices) 



m^y (2 + //)2 - z/ r ^(2 + /x)+zy n{n 



2 



= 87r(l + /.)3/2 + 2(1 + ^) ^3 + - ^m^e ^ , 



^ /C2 -I- n)'^ — y 

.0+ = £^0+ + 16^(1+^)5/2 - -B2 + B^ + 2S4 - //mA4 - m(2 + //)A6} 

(2.11) 



with the Ai{s,u) and Bi{s,u) evaluated at threshold sth = rn?{l + A*)^ and 
■J^t/i = "7,^(1 — 11 — iJ,'^ + A*^)/(l + A*)- This completes the necessary formalism. 

III. EFFECTIVE LAGRANGIAN 

In this section, we will briefly discuss the effective lagrangian on which our one-loop 
calculation in CHPT is based. We heavily borrow from our previous work on threshold 
pion photoproduction with = 0. Here, we will only discuss the terms which have to be 
added because of A;^ < in the electroproduction case. For the details we refer the reader 
to I. 

To systematically work out the consequences of the spontaneously broken chiral sym- 
metry at low energies, one makes use of an effective lagrangian of the asymptotically 
observed fields, in our case the Goldstone bosons (pions) and the nucleons {i.e. proton and 
neutron). We work in fiavor SU(2) and mostly in the isospin limit m„ = = m. The 
proton and neutron are collected in the isodoublet Dirac-spinor, 

*=(:) (3.1) 



which transforms non-hnearly under the chiral group SU(2)l x SU{2)r, 



^{x)^IC[gL,gR,U{x)]^{x) (3.2) 

where gL,R are group elements of SU{2)l^ii and the unitary unimodular matrix U{x) 
embodies the pion fields. For our calculation, it is most convenient to work in the so-called 
cr-model gauge, 

U{x) = [a{x)+iT-TT{x)]/F, a^ + jr^ = F^ (3.3) 

with F the pion decay constant in the chiral limit. To leading order, one calculates tree 
diagrams from the effective lagrangain with the least number of derivatives and quark mass 
insertions. It is given by 

4'i + 4^i (3.4) 

where the upper index indicates the so-called chiral power according to the counting rules 
spelled out in I. The explicit expressions for the lowest order pion-nucleon, /^^^j^, as well 

(2) 

as the pion lagrangian, Onr, are also given therein. At next-to- leading order, one has two 
types of contributions. These are the one pion loop diagrams which only involve the few 
parameters occurring at lowest order. However, there are also local contact terms which in 
general are necessary to renormalize the divergences of the one-loop graphs. To one loop 
order, i.e. chiral power q^, their generic form is 

4i + + A/:l°i + A£« + 41? . (3.5) 

The terms AC^^j^'^ relevant to renormalize the nucleon mass and pion-nuclcon coupling 
constant are discussed by GSS [11] and in I. In the meson sector, renormalizes the 
pion decay constant F-^ and the pion mass and there is one term which will give extra 
contributions due to < 0. It enters the pion charge form factor [5], 

4'i = ^|^(K,«1/;.) (3.6) 

with 

f^, = e{d^A,-d,A^){uQu^±u^Qu), 

(3.7) 



le = -T+ ^6(A), L = ^\ + _(^^ _ 1 _ in47r) \ (3.9) 



Here, Q = diag(l,0) is the charge matrix and denotes the (external) photon field. 
The constant Iq has a divergent piece (pole at d=4), its finite part can be fixed from the 
empirical value of the pion mean square charge radius, 

where the following renormalization prescription was used to cancel the infinities from the 
loops 

with A a scale introduced in dimensional regularization. Using the empirical value 
< >^= 0.439 fm^ [12], we have Zg(l GeV) = 6.6 • IQ-^. The other low-energy constants 
which occured already in photoproduction are I3 and I4. These are related to the chiral 
corrections of the pion decay constant and pion mass at order m and m^, respectively. In 
the ttN sector, there are three new terms contributing to pion electroproduction at order 
q^. These read 

+ ff ^>l3*757'^([i^^/^-J + ^K, 

The first two terms in eq.(3.10) are related to the electric mean square charge radii of the 
proton and the neutron. It is well-known that these radii develop a logarithmic singularity 
in the chiral limit. The diagrams responsible for this behaviour are in fact divergent and 
the constants 69, bg absorb these infinities according to the renormalization prescription 

= f (1 - 9\) + bliXl h = ^{9\-l) + m) . (3.11) 

Notice that the combination 69 + 269 is finite. The pertinent Dirac form factors of the 
proton and neutron to one loop order read, 

+ ^2^{-4%(fc2) - 16m2^4(fc2) - 2f 3(fc2) + 4m2f 4(fc2) + k^it^ie) - t^{e)]] , 



= 7777 . 



+ gl^{A^^{k^) + 16m2^4(A;2) _ A^ik^) + Sm^f 4(A;^) + 2k^[f4{k^) - f^{e)]]. 

(3.12) 
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where the loop functions 7i(t) and rj(t) are defined in appendix B. From the empirical 
charge mean square radii, < r% >p= 0.74 fm^ and < r% >„= —0.12 fm^ we can fix the 
coefficients at A = 1 GeV as 6^(lGeV) = 3.85-10-3 and 6^(lGeV) = 1.45-10-^. The last 
term in eq.(3.10) is related to the slope of the axial form factor of the nucleon, To 
one-loop order it can be written as 

GA{t) = ^^|l + ^&i3 + qIn [-2f 3(t) - 4m2f 4(t) + t{t4t) - t^it))] I (3.13) 

The one loop contribution to the axial mean square radius is finite and very small, therefore 
the low-energy constant 613 amounts for most of its value. Using the average empirical 
value from (anti)-neutrino scattering experiments, < >= (0.416 ± 0.02) fm^ [13] we 
find 613 = (10.05 ± 0.62) GeV-^. All other terms in jC^^^ which are relevant and already 
appeared in the photoproduction process are discussed in I. 

IV. INVARIANT AMPLITUDES TO ONE-LOOP 

In this section, we will be concerned with the chiral expansion of the invariant ampli- 
tudes of pion electroproduction as introduced in section II. To keep this section short we 
will give only the expressions for the reaction 7*^ — > n^p which was already discussed in 
the letter [14]. In appendix A, the amplitudes for 7*n n^n as well as the minus am- 
plitudes A[~\s,u),B^~\s, u) can be found. Together with the ones displayed here, these 
allow to calculate the charged electroproduction channels 7*p n^n and 7*72 'K~p 
by using eqs.(2.7). To keep the notation compact, we will give ^i,3,4,6('S, w), B2{s^u) and 
Bi{s,u) + 234(3, u) or an independent subset of the Bi{s,u), {i ^ 3,5). These are then 
converted into the Ai{s,u) via eqs.(2.6). 

IV. 1. GENERAL REMARKS 

We are seeking the chiral expansion of the invariant amplitudes Ai{s,u), i.e. the ex- 
pansion in small external momenta and quark masses, 

Ai = AY^' + Al-^'"'^ + Af . (4.1) 

Here, Af^^^ is the contribution calculated from the lowest order effective lagrangian eq.(3.4). 
^i-loop g^gjj^g fj-Qj2i the one pion loop diagrams generated by the lowest order vertices 
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and involves only the parameters gA, F and m and finally is calculated from tree 
diagrams with exactly one insertion from or C^^^ ■ They carry the information from 
the local counter terms whose coefficients are not fixed by chiral symmetry requirements. 
Before giving the explicit formulae, it is instructive to first discuss the general structure 
of the Ai{s, u) which emerges at one loop. As in the case of photoproduction there are 66 
topologically inequivalent one-loop diagrams. These can be grouped into three separately 
gauge invariant classes, labelled class I, II and III, respectively. While class I scales like qa ~ 
Qt^n (if one uses the Goldberger-Treiman relation), the diagrams in class II and III scale 
as g^j^i cf. figs. 2, 3, 4 in I, with Qt^n the strong pion nucleon coupling constant. In essence, 
any of the (s, u) has a contribution from electric and magnetic Born terms, the non-pole 
loop contributions from classes and III plus the additional counter term contributions. 
The electric Born term is the tree amplitude multiplied with the corresponding hadronic 
electromagnetic form factor consistently calculated to one-loop and it involves the physical 
values of the pion-nucleon coupling Qt^n and the nucleon mass m. For the channel 7*p — > 
7r°p, the relevant form factor is Ff{k'^), the Dirac form factor of the proton. Similarly, 
the magnetic Born term involves the proton Pauli form factor F2 (A;^) as generated by 
some graphs in class II and a counter term which allows to adjust the proton anomalous 
magnetic moment. Notice that the latter is not present at tree level but rather generated 
by loops and higher order counterterms. As in the photoproduction case many of the one 
loop diagrams simply lead to a mass and coupling constant renormalization (as discussed 
in I). Another important topic concerns the isovector nucleon electric form factor F^{k^) 
and the pion electric form factor F^{k^). These appear naturally in some of the amplitudes 
exhibited below. Since we are starting from a gauge invariant theory and because these form 
factors are build up by loop and counterterm contributions in a fashion that conserves the 
pertinent symmetries, there is no need of setting F^^k"^) = F^{k'^) as it is done in many 
model calculations to preserve gauge invariance [2]. This feature is particularly clearly 
shown in the discussion of the low-energy theorems in section 5. After these general 
remarks, let us now discuss the various contributions to Ai(s, u) for the process 7*p — > 7r°p 
at one loop order. 

IV.2. ELECTRIC AND MAGNETIC BORN TERMS 
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The electric Born term modified by the proton Dirac form factor F[{k^) (see eq.(3.12)) 
leads to the following amplitudes 

Ai{s,u) = -B2{s,u) = eg^Nl — ^ + ^ ) Fl{k^) , 

\s — u — J (4.2) 

Si(s, U) + 2^4(5, U) = ^3,4,6(5, W) = . 

It is important to notice that Ff{k^) is not the full physical form factor but rather its one 
loop expansion. This is a point frequently overlooked in the literature. We will come back 
to this point in section V in connection with the discussion of the low energy theorems. A 
brief discussion of this point has already been given in ref.[14]. Part of the magnetic Born 
term comes from the counter term which corrects the proton anomalous magnetic moment, 

^'^^'''''"2m2 

= - '-^^"^ ( - ] ' (4-3) 



'P1 


^2,5,6' 




' 1 




— m? 
' 1 


6Kp\^ 




, s — m? 



where 



A^{s, u) =- -^oKp { - — ^ + 



5Kp = ^vp + 8^^^m2{^4(m^0)- 274(m^0)} = 1.288 (4.4) 



is adjusted to give the empirical value of the proton anomalous magnetic moment, Kp = 
1.793. In contrast to what was done in I, we have not subtracted the loop contribution to 
F2 (A;^) from class II and therefore only Sup appears in eq.(4.3). 

VI.3. CONTRIBUTIONS FROM CLASSES I, II, III AND COUNTER 
TERMS 

For class I and II the following representation is most economic 

Ai{s,u) = ai{s, k"^) + ai(w, /c^) - B2{s,u), 

B2is, u) = — — 2 [biAim'^, k^) - &i4(s, k"^)] H — 2 [^i4("^^; k^) - bi4{u, k^)] , 

Bi{s, u) + 2i?4(s, u) = 614(5, A;^) - 6i4(w., k^)-, ^■^{s, u) = a^i-s, k'^) - a^iu, k^), 
A4^{s,u) = 03(5, /c^) + a3{u,k'^), A6{s,u) = ae{s,k'^) - a^iu^k^) . 

(4.5) 
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Class I, which scales as Qt^n-, leads to the following contributions, 
ai(s,/c^) = ^^J^(s-m^)[76(s,/c^) -74(s,A;^)], 

TT 

614(5, k^) = ^{(5 + 3m2)[276(5, k^) - 72(5, A:^)] + {s - m')[7i(5, k") - 275(5, A:^)]}, 

TT 

03(5, /e^) = ^^1^77174 (s,/c^), 

TT 

a6(s, /c^) = ?^J^m[276(s, /c^) - 72(5, k"^)] . 

(4.6) 

The various loop functions 7i(s, /c^) are defined in appendix B. 

The class II diagrams give the following expressions proportional to g^^, 

ai{s, e) = ^{{s- m2)[474(s, /c^) - 475(5, k^) + T^{s, k^) - T^is, k^) 

-2T4{s,k^) +2T6{s,k'^)]+8m^[2'y4{s,k'^) -T4{s,k^)]}, 

3 

614(5, = f^{(5 + 3m2)[475(5, e) - 271 (s, + ri(s, _ r2(5, /e^) _ 2r5(5, k^)] 



4m2 

+ [272(5, k') - 476(5, k') + 2r6(5, /C^)]}, 

s — 

3 

03(5, = !|EiV |4^g(,^ ;.2) _ 8^^(^^ ^2) ^ 4P^^^^ ^2) _ p^^^^ ^2) ^ p^^^^ ^2) 

- 2r6(5, + ^^[r4(5, k^) - 274(5, k^)]}, 
s — m^ 

3 

ae{s, e) = ^{475(5, _ 271(5, k'') + Ti{s, k") - T2{s, k") - 2T^{s, k") 



2m 



^^^±^[72(5, e) - 276(5, e) + r6(5, e)]] . 



+ 

s — 

(4.7) 

Notice that the one loop representation contains the magnetic (Pauli) form factor of the 
proton through the k"^ dependence of the loop functions. It is proportional to the residue 
of 03(5, k"^) at s = m^. Furthermore, one observes that 05(5, k"^) and 614(5, k"^) do not have 
such poles since 27fi(m^, fc^) — 72(m^, /c^) = = r6(m^, fc^). This can be easily seen from 
the Feynman representation eq.(B.2,3) by changing the integration variable x ^ 1 — x. 
Finally, we have to give the amplitudes for class III, which are most tedious to work out. 
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A compact representation is obtained by giving the Bi{s,u) {i ^ 3,5), 



M2 

n 

s — 

J2 n^2\ 



B,{s,u) =egl^{ rj^(5) - —^r-^{s) - rj^(s) + {s-m')G^{s,t,k') - 2G^{s,t,k') 



- k^G5{s, t, r) + {u- s)G6{s, t, r) + (t + r - m^)G7{s, t, r) - tG^is, t, r) 

2 

+ {t- 4:m^)Gs{s, t, + 2n^{s, u, k'') + [474(5, P) - 475(5, P) - Ti{s, k") 

+ T2{s, e) + 2T^{s, e) - 2Tq{s, P)] + ^ [-271(5, k^) + 272(5, P) + 475(5, P) 

- 476(5, k") - ri(5, k") + r2(5. A;') + 2r5(5, k") - 2Vg{s, k")] | + (5 ^ , 

i?2(5, =e^^^|-r-^(5) + — ^£-^(5) + r-^(5) + t,(t) + (5 - m") [G2(^, fc') 

- ^3(5, t, /c^) + 1^2(5, u, k^) - 2ne{s, u, k"^)] + 2G^{s, t, k^) + k^Gr,{s, t, k^) 
+ (2m2 - 5 - u)Gq{s, t, k^) + {M^-t- k'^)G7{s, t, k'^) + ^^9(5, t, k^) 
+ {u- s-Am^ + Ml- k^)G8{s, t, k^) + 2(s - - M^) [^10(5, k^) 
+ flio(5, /c^)] - 2114(5, u, k^) + 2{Ml - 3m^ - 5)118(5, /c^) 
+ 272(5, + _^ [27i(s, _ 272(5, e) - 475(5, fc2) + 476(5, e) 

+ ri(5, k^) - r2(5, A;2) - 2r5(5, k^) + 2r6(5, A;^)] | + (5 ^ «) , 

5i(5, «)+254(5, «) = e5r^^|(5 - rr?) [Gi(5, t, k^) + ^2(5, t, A;^) - ^3(5, t, A;^) - 2^5(5, t, A;^) 

+ Oi(s,ii, A;^) -03(s,w,A;2) - 2^5(5, w, A;^)] + 2(M2 - Sm^ - s) [^0(5, t, A;^) 

- Gio(5, t, A;^) + 06(5, -u, A:^) - Oio(5, u, k^)] + 2(25 - 2m2 - M^) [07(5, t, k^) 
+ 07(5,'u,A;^)] +2(m^ + M^ -5)[G9(5,t,A;2) + 09(5,'u,A;^)] +7o(5,A;2) 

+ 272(5, e) + 475(5, e) - 476(5, A;2)] I - (s ^ «) , 

^6(5, u) = -^7^^{ri(5, A:2) - r2(5, e) + 2r6(5, e) + 475(5, A;2) 

- 8m^G8(s, ^, A;^) - 8m^n8(s, A;^)} + (5 ^ it) , 
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Bris, u) = ^^3^{ri(s, e) - r2(s, k^) - 2r5(s, k^) + 271(5, e) 

— 475(5, /e^) - 8m2G'6(s, /c^) + 8m2Gio(s, /c^) 

— 8m^ilQ{s, u, k"^) + 8m^f2io(s, tt, /c^)} — (s -^-^ tt) , 

58(s, u) = ^^3^{ri(5, A;2) - V2{s, e) + 2r6(5, A;2) + 476(5, fc^) 

— 8m^Gio(s, t, /e^) — 8m^riio(s, it, /c^)} — (s •(-^ it) . 

An excellent numerical check on these rather involved expressions is given by the gauge 
invariance relations eq.(2.3) (for that, one also calculates S3 5(5, w)), with the exception 
of Bi{s,u) which is not constrained by the gauge invariance condition k^J^ = 0. As a 
further check, all the above expressions of course match to the formulae given in I for the 
case k!^ = 0. Finally, for the channel 7*p n^p there is one counter term contribution 
to ^4(5, u) — aj from C'j^}^ The same term was already present in the photoproduction 
calculation and its coefficient was estimated from resonance saturation. Equivalently, one 
can say that the determination of the constant ^ in I amounted to an overall best fit to 
the total photoproduction cross section for — > ir^p in the threshold region based on the 
Mainz data [15] (with < 160 MeV). 

With the invariant amplitudes given here and in appendix A, we are now in the position 
of discussing the pertinent physics issues related to threshold pion electroproduction. 

V. LOW ENERGY THEOREMS 

V.l. GENERAL REMARKS 

Pion electroproduction low energy theorems (LETs) have been discussed in many 
articles [16]. One of the most recent discussions is due to Scherer and Koch [17]. The 
underlying idea is to expand the threshold S-wave multipoles £'0+ and Lo+* in terms of 
the dimensionless and small parameters fi = MT^/m ~ 1/7 and v = k'^/m?. Evidently, v 
can be made arbitrary small as k'^ approaches zero, the so-called photon point. One of the 
main differences to the photoproduction case was first pointed out by Nambu, Lurie and 

* Here £'0+ and Lo+ stand as generic symbol for any isospin channel. 
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Shrauner (NLS) [18], namely that charged electroproduction involves the axial form factor 
of the nucleon Ga (t) ■ It therefore gives an other possibility of determining this fundamental 
quantity. However, the LET of NLS refers strictly to the chiral limit and one therefore 
has to add corrections for the physical case of a finite pion mass. Such corrections are 
discussed for example in the monograph [2]. As we will show in some detail later, CHPT 
allows to systematically calculate the next-to- leading order corrections to the LET of NLS, 
as already reported in ref.[19]. This is only one particular example of the chiral expansion 
done consistently. In fact, the one loop calculation performed here allows us to calculate 
all terms of the multipole amplitudes £'o+(/i, i') and Lo+(/i, ly) up to and including order 
0{ii^, u). The results of this calculation will be reported below. 

First, however, some clary fying words about the meaning of LETs and their relation to 
the CHPT calculation are in order since in the present literature one finds many erroneous 
statements. In QCD, the quark masses are a priori free parameters. For the two fiavor 
sector, one can set 111^ = rrid = to a good approximation. In this limit, the so-called 
chiral limit, the pions are exactly massless, as mandated by the Goldstone theorem. CHPT 
now allows to systematically work out the consequences of the spontaneous and the explicit 
chiral symmetry breaking in QCD, as an expansion in small external momenta and quark 
masses. For the case of pion electroproduction at threshold, these small expansion param- 
eters are the pion mass ~ m and the photon four-momentum k^. As it is well-known 
and most clearly spelled out by Weinberg [20] , CHPT embodies the very general principles 
of gauge invariance, analyticity, crossing and PCAC (pion pole dominance). Therefore, to 
lowest order, one recovers the venerable current algebra LETs, which are based on these 
principles. The effective chiral lagrangian is simply a tool to calculate these but also all 
next-to-leading order corrections in a systematic and controlled fashion. The question of 
interest here is now: What are the corrections to the lowest order statements for £'o+(//, v) 
and Lo+(//,z^) at threshold? This question can be unambigosly answered by making use 
of the CHPT machinery. For that it is mandatory to consider pion loop diagrams. It was 
already shown in the photoproduction case that the conventional LET due to Vainsthein 
and Zakharov [21] and de Baenst [22] was incomplete and has to be modified at next-to- 
leading order, 0{p'^), due to a logarithmic singularity of certain loop diagrams in the chiral 
limit of vanishing pion mass. Such an effect can never be found without explicit calculation 
of all relevant diagrams. In the light of this result we also expect similar modifications to 
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show up in pion electroproduction. This expectation was aheady borne out by an exphcit 
calculation of the class I diagrams for the reaction — > 7r°p [14] . In the following section, 
we will give the complete list of LETs for electropion production (in the isospin basis). 

In view of the above arguments, it should be clear that simply calculating Born (tree) 
diagrams supplemented by form factors (in a gauge invariant fashion) can not give all 
corrections at next-to-leading order [17,23,24]. To be more specific, consider first charged 
pion electroproduction. The leading term in £'o+(/i, i^) and Lo+(/i,^') are of order one, 
which is nothing but the generalization of the famous KroU-Ruderman theorem [25]. In 
ref.[17] it was claimed that all corrections up to and including v) can be calculated 

from PCAC and current conservation without considering loop diagrams. From the above 
remarks on a consistent expansion in CHPT, it should be obvious that this expectation 
is too naive. It is based on the incorrect assumption that while PCAC is a very general 
principle, CHPT is derived from a particular lagrangian and would be more restrictive. 
Quite contrary, any theory or model which claims to embody PCAC should lead to the 
same result as CHPT if one goes to the same level of sophistication. Similarly, the recent 
discussion of Ohta [24] is misleading. When calculating the LETs in QCD, one never 
has to worry about half ofi-shell nucleon form factors and thus the findings reported in 
refs. [14,19] are model-independent and do not rely on any assumptions about these form 
factors. We hope that with these remarks we have been able to shed some light on the 
seemingly controversal theoretical interpretation of the LETs. In fact, the meaning of the 
LETs is unique (experimental quantities are expanded in the chiral symmetry breaking 
parameter M^r and small external momenta) and their accuracy at a given order can be 
tested experimentally. Let us now present the pertinent complete results at next-to-leading 
order. 

V.2. FORMULATION AND DISCUSSION OF THE LOW ENERGY THEO- 
REMS 

The derivation of the LETs for E^'^'~^ and Lq"^'*^'"^ in electropion production is a 
straightforward generalization of the photoproduction case investigated in I. Here, we only 
wish to point out that for this calculation it is most convenient to make use of the heavy 
mass formulation of baryon CHPT [26]. It was shown in ref.[27] that a one loop calculation 
in this framework reproduces all relevant terms. We therefore skip the computational 
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details here. The correct form of the LETs for Eo+{iJ,,u) and Lo+{iJ,,u) in the various 
isospin channels (+, 0, — ) at next-to-leading order is given by* 



-(+) 

"0+ 



-(-) 

"0+ 



327rm 



327rm 



{-2;u + ;u2(3 + 



327rm 



87rm 



87rm 



V(2 + /^)^-^ 



2(l + /x)3/2(jy-2/x2_^3) 



(5J) 

where the functions E]j(— z//|U^), (j = 1,2,3,4) can not be further expanded since the 
argument — counts as order one. They are given by 



2i(p) = 

S2(P) = 
S3(P) = 
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+ TTT^ r-;77^ arccos 
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S4(p) 



+ 2 / (ia; \/(l — a;)[l + a;(l + p)] arctan — == 
Jo x/(T^ 

— / da; — , 



a; 



V(l-a;)[l + a;(l + p)] 



a:(l - 2a:) 



arctan 



X 



V(l-a;)[l + a;(l + p)] ^ (i - + x{l + p)] ' 



These functions are shown in fig.[l] for < p < 10. They exhibit a very smooth behaviour. 

3.71 and Kg = Kp + Kn = —0.12 are the isovector and isoscalar anomalous 
magnetic moment of the nucleon, respectively. We have to make various remarks about 



* Some of the kinematical prefactors in eqs.(5.1) have not been expanded to keep the 
notation compact. 
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the LETs exhibited in eqs.(5.1). First, adding the (+) and (0) multipoles, one recovers 
the LETs for 7*p n^p derived in refs. [14,27]. One sees that class II and III do not at 
all contribute at order 0{ji'^^v). This feature is most easily understood in the heavy mass 
formulation of CHPT, where simple selection rules determine those few diagrams which are 
non- vanishing at threshold (see ref.[27]). Second, notice that the LETs do not contain the 
full electromagnetic form factors of the nucleon. This is due to the fact that in the consistent 
chiral power counting on which CHPT is based these form factor effects are already of order 
0{jiv) = 0{q^) and therefore consistently have to be dropped. Furthermore, one sees only 
the normalization of the magnetic form factors, i.e. the respective anomalous magnetic 
moments. Of particular interest are the LETs for the (— ) amplitudes. In the chiral limit, 
the LET for Eq^^ agrees with the one of NLS [18] at order 0{u). The constant C which 
appears also in photoproduction is of order one and it is given by 

C = ^5^(l + ln.) + ^{4-^-'+4-'} (5.3) 

where the coefficients a[~^ have been estimated in ref.[l] using resonance saturation. For the 
central values given there, a[^'' = 1.4 GeV~^ and a^~^ = -9.9 GeV~^, we find C = 0.40. In 
ref.[19], the consequences of the LET for E^'^^ were discussed. If one expands S3(— 
in powers of = vw? and picks up all terms proportional to v [i.e. all terms which 
contribute to the slope of E^^ at A;^ = 0), one finds 



dE, 



(-) 

'0+ 



where the first two terms on the right hand side of eq.(5.4) have first been given by NLS. 
As discussed in ref. [19], the new term at order k"^ (not vanishing in the chiral limit) due 
to loop effects has a numerical value of —0.0456 fm^ and this closes the gap between the 
experimental determination of the nucleon axial mean square radius < t\> from neutrino 
experiments, < r\ >^= 0.42 fm^, and pion electroproduction analysis, < r\ = 0.37 
fm^ [31]. To get an idea about the higher orders in k^, we show in fig. [2] the two functions 



$i(fc^) = l + y <ri>. 



^ n2^ 2 \^ I k\ 7r2 1) 

^.(k^) = l + -<rl> +^|-3(-^) - - - -| 



(5.5) 
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for < r\ >= 0.42 fm^. For < 0.5 GeV^, the corrections of order /c^ (and higher) from 
one loop are quite smaU. Most important for the analysis of new precise electroproduction 
experiments is that one takes into account the novel term at order appearing in eq.(5.4). 
Let us now turn to the LET for Lq~^ — Eq^ in eq.(5.1). In contrast to the electroweak 
nucleon form factors, the full one loop expression of the pion charge form factor F^{k^) 
enters. This can be traced back to the following fact. The one loop representation of 
F^{P) takes the form 

F^ik') - 1 = -^jci/c^ InM. + c.^ + + . . .| (5.6) 

and all terms in the curly bracket are of the same order 0{q'^). This feature is similar to 
the functions E]j(— 1^//~^) whose expansion can not be truncated at any finite order. To 
the order we are working, one can even identify the one loop pion charge form factor with 
the empirical one, since the differences are order g^. To get an idea of the higher order 
corrections, we can use the two loop representation of ref.[28]. At k^ = —0.2 GeV^, i.e. v = 
—0.23, one finds that the one loop and the two loop representation of F^{k'^) differ by 
22%. Finally, let us stress again the two salient features of the LETs presented in eqs.(5.1). 
First, in some cases virtual pion loops generate contributions with chiral singularities which 
modify the form of the LETs based on an incomplete calculation of tree diagrams including 
electroweak form factors. These unfamiliar terms are given by the S-functions. Second, in 
the consistent expansion only the first moment of the nucleon electromagnetic form factors 
and the first and second moment of the nucleon axial form factor survive. The pion charge 
form factor, however, enters with its full one loop expression. The LETs presented in 
eqs.(5.1) are the ones which follow from the spontaneous (and explicit) beaking of chiral 
symmetry in QCD. 

It might be instructive to investigate how fast the convergence of the LETs is. In 
the photoproduction case it was shown that the 0{jJL^) contributions of the one loop ap- 
proximation to £'o_|_^(/i) substantially reduce the 0{iJ?) terms. Such a trend can also be 
seen here. To be more specific, consider Lq^{ii,Q) at the photon point. For that we have 
evaluated all the class I contributions up-to-and-including 0{ijfi, In//). The result reads 



" " ' ' " -j.-r — UL ill LA, —r ub \ — — — — — I z' 

(5.7) 



T-^ V( n\ '^yvrA' I 1^2 -^^ 3 I 
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= -0.0247 • (1 - 0.797 + 0.377) GeV"^ = -0.0143 GeV"^ 



20 



where we have exhibited the resulting contributions at order (U, //^ and fj,^, in order. One sees 
that the series actually converges slowly for fj, ~ 0.14. Notice, however, that these results 
should be considered indicative since higher loop corrections will modify the coefficients 
of the fjfi and //^In/x terms. The result of the full one loop calculation will be given in 
section VI and compared to the recent experimental determination [3]. In the charged 
channels the corrections are increasing with with increasing \k^\ due to form factor effects 
(see section 6). This is different from the photoproduction case, where the loop corrections 
to E'^^~^^^'^{fi) and E'^^~''^ ^{jj) are small and move the chiral prediction closer to the 
empirical values. 

V.3. INCLUSION OF SOME ISOSPIN BREAKING EFFECTS 

As already stated, the LETs presented so far refer strictly to the isospin limit and do 
not account for isospin breaking either through quark mass differences (~ rud — my) or 
(higher order) electromagnetism (~ e^). At present, we are not able to implement in a 
consistent way all possible effects of isospin breaking. However, to get an idea about the 
size of such effects, we will present the results of a simplified approach. It is well-known 
that the mass of neutral and charged pion differ slightly. 



0.033 (5.8) 



and this difference is almost entirely of electromagnetic origin [29]. One therefore expects 
some effect to come from this mass difference of pions in the loop and the asymptotic 
state. Indeed, taking this mass difference into account in neutral pion photoproduction 
leads to a substantial correction and Eq^~~''^ ^ = —1.97 • 10~^/M^+ in good agreement 
with the commonly accepted empirical value Eq^"''^ ^ = —(2.0 ± 0.1) • 10~^ / [7].* 
A similar mechanism is also operative in the function Si 2(p) which enter neutral pion 
electroproduction. Expanding in hq = M^o/m and using po = —^/l^'o we find 

Si(Po) = ^|-\/m2_^-M,%+£ dx^M^^-x^Mi„ + ex(x-l)\ , 

° - ^^^^0 + h'^xix - 1) 



* In ref.[7], isospin breaking was only considered in the class I diagrams which give the 
dominant loop-effect at /c^ = 0. 
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In fig. [3], we show the functions Si^2(po)- Their dependence is similar to the isospin 
symmetric ones =1^2 (p), but the value at the photon point is substantially reduced. While 
when neglecting the tt^ — 7r° mass difference, Si(0) = tt ~ 3.14 and S2(0) = 2 — tt ~ —1.14, 
we have now Si(0) = 2.28 and S2(0) = —0.74, which amount to a sizable reduction 
of 27% and 35%, respectively. The main point to stress here is that in neutral pion 
electroproduction it is mandatory to go beyond 0{q^) and to include at least in a simplified 
manner isospin breaking effects (like the charged and neutral pion mass difference). A 
systematic study of these effects is beyond the scope of this paper. 

Differentiating between the neutral and charged pion masses in the loops amounts to 
a redefinition of some loop functions as detailed in appendix D. In essence, one has to 
disentangle the two thresholds, s^^^^ = m + M^^o and s^^. = m + Mt^+ ( the neutron-proton 
mass difference is not accounted for). In particular, some of the low-cncrgy constants which 
enter the proton charge radius and anomalous magnetic moment have to be readjusted. We 
find 6g(lGeV) +69(lGeV) = 5.17-10""^ as compared to 5.30 ■ 10""^ in the isosymmetric case 
(see eq.(3.12)). Similarly, 5Kp defined in eq.(4.4) changes from 1.288 to 1.291. That this 
prescription conserves gauge invariance can be checked after constructing the amplitudes 
Sj(s, u) and using eqs.(2.3). We should also stress that in this case of pion electroproduction 
one has to include the pion mass difference in the diagrams of all classes (the procedure 
used in ref.[7] for photoproduction is, however, legitimate). 

VI. RESULTS AND DISCUSSION 
VI.l. ISOSYMMETRIC CASE 

First, we consider the isosymmetric case, i.e. we do not differentiate between the 
proton and the neutron masses and also not between the neutral and charged pion masses. 
Throughout, we set rup = rUn = m = 938.27 MeV, = 93.1 MeV, fir^jv/^TT = 14.28 
and e^/47r = 1/137.036. When considering charged pion production we set = M^± = 
139.57 MeV and for neutral pion production = M^o = 134.97 MeV. This is a consistent 
procedure within the one-loop approximation we are using. 

In Fig. 4 we compare the CHPT prediction with the recent NIKHEF data [3]. In the 
same figure, we also show the results from the tree diagrams (open circles) and from the 
pseudovector (PV) Born terms with form factors (open squares) calculated in Ref.[8]. Note 
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that all of the PV results presented here do not include vector meson-exchange and final 
state interaction terms introduced in [8] . We see that the one-loop contribution drastically 
reduces the calculated cross sections and brings the CHPT results within the experimental 
uncertainties. The results from the PV terms are also not too different from the data. To 
have a more rigorous test of CHPT more precise measurements at small k"^ are needed. 

Let us now turn to analyzing the dynamical content of the CHPT calculations. We 
consider a kinematical situation with the final pion-nucleon invariant mass W = 1080/1074 
MeV (for charged/neutral pion production) which is perhaps most realistic for future exper- 
iments. At this near threshold energy (Wthr = 1077.84/1073.24 MeV for charged/neutral 
production) the cross sections are dominated by the S-wave amplitudes £^o+ ^ind Lo+ 
(5'o+)- In Figs. 5 and 6 we display the real parts of these two amplitudes for p(7*,7r+)n, 
n(7*, 7r~)p and p(7*, i^^)p, in order. In all cases the tree and PV results arc rather similar 
whereas the one-loop prediction is significantly different from these. This comes in part 
from the irreducible class I, II, III diagrams and to a large extent from the difference in the 
/c^ dependence of the pion form factor and the isovector nucleon charge form factor which 
is a natural ingredient in our approach.* In Ref.[30] we had already shown that the S-wave 
cross section for neutral pion production off the proton 

ao = l-Eo+P - e-2 |Lo+P (6.1) 

where e and = {s — rn^ + k^^jl^ represent, respectively, a measure of the transverse 
linear polarization and the energy of the virtual photon in the ttA^ rest frame, follows 
nicely the data of Welch et al. [3] whereas PV and tree are at variance with the data with 
increasing In this channel we find at the photon point (/c^ = 0) |-^o+P = 0.2 fib in fair 
agreement with the result of [3], l-C'o+P = 0.13 ± 0.05 Clearly, a precise measurement 
of the k'^ dependence is the most effective way to test the CHPT predictions. We should 
stress that one should consider \k'^\ < 0.05 GeV^ since otherwise the loop corrections are 
too large. To further illustrate this point, we show in Fig. 7 the transverse and longitudinal 
cross sections ut and ctl, respectively, for p{'j* ,7r'^)n, n{'j*,7r~)p and p(7*,7r°)p, in order. 
In all cases the one-loop predictions are rather different from the PV and tree results. 
We find an enhancement of (7^ and a depletion of cry. Note, however, that the sums 

* The role of the various form factors has also been stressed by Olsson et al. in ref.[31]. 
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cr = ctt+cl are very similar for tree, PV and one-loop. To test the dynamical content of the 
chiral symmetry prediction, it is therefore mandatory to perform a transverse/longitudinal 
separation. To further exhibit these differences, we plot in Figs. 8 to 11 the differential cross 
sections daT/dfl, daL/dQ, daj/dfl and dap/dft for three values of k'^ = —0.001, —0.04 and 
—0.08 GeV^. The most striking difference appears in dar/dO, for neutral pion production. 
While the CHPT prediction becomes forward peaked for > 0.04 GeV^, the PV and tree 
results remain backward peaked until ~ 0.07 GeV^. For < 0.05 GeV^ the shapes 
of dai/dfl, and dap/dO, are not too different between the tree, loop and PV predictions 
in case of charged pion electroproduction. For neutral electropionproduction, the angular 
distribution dap/dfl shows a significant difference between the tree, PV and one loop 
results even at low photon four-momenta. We will come back to this in section 6.2. It is 
worth to stress that as increases, the contributions from the diagrams of classes II and 
III become more and more important. To illustrate this, consider the S-wave cross section 
(6.1) at P = -0.06 GeV^. The reduction of the one loop result compared to the tree 
(or PV) prediction [30] is almost entirely due to the class II and III diagrams (with equal 
share). The class I diagrams contribute insignificantly to ao at this k^. This is completely 
different from the photoproduction case discussed in refs.[l,7]. 

It is important to note that these results at W — 1074/1080 MeV shown in Figs. 5 to 
11 contain contributions from P-wave multipoles. In the threshold region these are small 
but they can nevertheless contribute significantly to the various cross sections through the 
interference with the S-wave amplitudes. In table la,b,c we give the real and imaginary 
parts for Eo+, So+, Ei+, Si+, Mi+ and Mi_ for P = -0.001, -0.04 and -0.06 GeV^ 
for p(7*, 7r"'')n, n(7*,7r~)p and p(7*,7r°)p, in order. The knowledge of these multipoles 
allows one to construct any observable one wishes to measure, in particular the many 
which appear in the case of polarized electrons. Some pertinent formulae are summarized 
in appendix C. 

To explore the sensitivity of the charged pion production amplitudes to the axial form 
factor, we have performed a set of calculations at fc^ = — 0.04 GeV^ varying the dipole mass 
between 0.96 and 1.16 GeV, corresponding to squared axial radii between 0.35 and 0.51 fm^ 
(remember that our central value is Ma = 1-06 GeV). From the LETs discussed before, 
eq.(5.1), one expects that a smaller axial radius leads to an increase in transverse strength 
and a decrease of the longitudinal strength. This expectation is borne out by the actual 
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calculations. In fig. 12 we show the differential cross sections dar/dQ, daL/dQ, dai/dQ and 
dap/dfl for this range of axial radii for p{'y*,7r~^)n (the results for 77,(7*, 7r~)p ^^e similar 
and not shown). Particularly sensitive to Ma is the ratio ar/cTL which decreases/increases 
by approximately 20 per cent for Ma lowered/enhanced by 0.1 GeV. In contrast, the sum 
cl+ct is almost independent of the axial radius which again points towards the importance 
of a transverse/longitudinal separation. A summary of these results is given in table 2. 

Furthermore, we have investigated the influence of the counterterm di which already 
appears in the photoproduction reaction p{^,'!T'^)p. Reducing its strength by a factor 2, 
the transverse cross section is diminished since the electric dipole amplitude is smaller 
(at the photon point). For finite fc^, however, this trend reverses and eventually leads to 
enhancement of ax (by 50 per cent at /c^ = —0.04 GeV^). The longitudinal cross section is 
essentially unaffected by this counter term. The charged electroproduction amplitudes are 
insensitive to this finite contact term. We should stress again that its central value could 
be considered as coming from a best fit of the precise total photoproduction cross section 
data [15]. We will come back to this when discussing the inclusion of some isospin-breaking 
effects. We did not vary the two counter terms appearing in the charged photoproduction 
channels since their influence on the pertinent observables is fairly small. 

VI.2. INCLUDING ISOSPIN-BREAKING 

As discussed in section 5.3, for the case of neutral pion elelectroproduction it is manda- 
tory to consider the dominant isospin breaking effect which stems from the pion mass differ- 
ence. In what follows, we consider — > ir^p using = 139.57 MeV and M^o = 134.97 
MeV for the pertinent loop functions (see appendix D). In all kinematical factors the neu- 
tral pion mass enters. Before presenting the results, it is worth to stress that the gauge 
invariance conditions (2.3) are fulfilled within machine accuracy. 

In fig. 13, the S-wave cross section oq defined in eq.(6.1) is shown in comparison to 
the isosymmetric CHPT [30] and the PV result (for e = 0.58). Although the pion mass 
difference has a pronounced effect on Eq^{ij,,0) and Lo+(//,0) (e.g. |ivo+P = 0.17//6), at 
finite k"^ the two CHPT curves do not differ significantly. The conclusion of ref.[30] that 
loop effects are needed to explain the trend of the NIKHEF data is not invalidated by the 
inclusion of isospin-breaking. This is due to the dominating effect of the class II and class 
III diagrams at finite /c^. 
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To get a better handle on the possible isospin-breaking effects, we show in fig. 14 the 
angular distributions daL,T,p,i/ dO. at k"^ = —0.04 GeV^. Only in the case of dap/dO, one 
finds a significant change from the isosymmetric case. The height in the peak differs by 
almost a factor of three. In all other cases, the differences are marginal. This can also be 
seen from the comparison of ut and gl shown in fig. 15 for —0.1 < A;^ < GeV^. The 
largest sensitivity of these results stems from the low-energy constant di. Reducing its 
strength by a factor of two, one finds an enhancement of dar/dO. in forward direction (at 
^ = 0, the enhancement factor is about 1.7) and of dap/dVL (the peak height is increased by 
about 1.4). The other two angular distributions remain almost unaltered. In the S-wave 
cross section ao this amounts to an increase of about 21 percent. We conclude that the 
largest uncertainty of our predictions stems from the knowledge of the low-energy constant 
d^. 



ACKNOWLEDGEMENTS 

Wc arc grateful to Nimai Mukhopadhyay for some useful comments. This work is 
partially supported by the U.S. Department of Energy, Nuclear Physics Division, under 
contract W-31-109-ENG-38. 



26 



APPENDIX A: INVARIANT AMPLITUDES OF PION ELECTROPRODUC- 
TION 

Here, we will complete the list of invariant amplitudes for pion electroproduction by 
giving those for the reaction 7*n n^n as well as the (-) amplitudes. 

A.l. THE CHANNEL 7*n ^ 7r°n 

Electric and magnetic Born terms: 



Ai{s,u) = -B2{s,u) = -eg^Ni — ^ + ^ ) F^i^) , 



{A.l) 



s — w? u — w? 

Bi{s, U) + 2S4(s, U) = ^3,4,6(5, w) = . 

with the neutron Dirac form factor Fi{k^) given in eq.(3.12). 

-44(S, U) =— dKn ^ H ^ 

2m \s — u — 



where 



is adjusted to reproduce the empirical value of the neutron anomalous magnetic moment, 

Kn = -1.913. 

For class I and II the following representation is again most economic 
Ai(s,u) = ai(s, /c^) + ai(tt, A;^) - B2{s,u), 

B2{s, u) = — ^ [6M(m2, e) - bi^s, e)] + — ^ [6i4(m2, k') - bi^{u, k^)] , 

s — u — 

Bi{s, u) + 2i?4(s, u) = 614(5, k^) - bi4{u, /c^), As{s, u) = a^is, k^) - a^iu, /c^), 
A4{s,u) = 03(5, /c^) + 03(16, /c^), Aq{s,u) = aQ{s,k^) - aQ{u,k^) . 

{AA) 
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Class I: 



ai{s, = ^{s- m')[2^e{s, - 274(5, + T,{s, k^) - V^is, - 2V^{s, 
+ 2r6(s,A;2)], 

614(5, = ^{(« + 3m2)[276(5, _ ^2(5, + 2r6(5, e)] 

TV 

+ {s- m2)[7i(s, /c2) _ 275(5, e) + ri(s, _ r2(s, _ 2r5(s, A;^)] } 
03(5, A;^) = ^^J^m[74(s, k'^) + T4{s, k"^)], 

■K 

ae{s, e) = ^^m[276(5, /c^) - 7^(5, A;^) + 2V^{s, k^)] . 

TT 



Class II: 



(^.5) 



3 

ai(5, = ^{(s - m2)[274(s, /c^) - 275(5, /c^) - V^{s, k^) + r2(s, A;^) 
+ 2r4(s, A;^) - 2r6(s, A;^)] + ^m''[^^{s, k'^) + r4(s, A;^)]}, 

3 

6i4(s, A:^) = ^{{s + 3m2)[275(5, A:^) - 71(5, e) - T^{s, k^) + T2{s, A;^) + 2T^{s, k^)] 

+ 2 ^2(5, A;2) - 276(s, k^) - 2r6(s, A;^)]}, 

asis, k") = ^{276(s, k^) - 474(5, A;') - 4r4(5, A;') + ri(5. A;') - r2(5. A;') 
m 

+ 2r6(5, k^) - -^^[U{s, k") + 74(5, A;2)]}, 

3 

06(5, A;2) = ^{275(5, k") - 71(5, k") - ri(5, A;2) + r2(5, k") + 2r5(5, A;2) 
m 

+ ^ii±^ [72(5, - 276(«, k^) - 2r6(5, e)] } . 

S — 777/ 



{A.Q) 
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Class III: 



We give the Bi(s, u) with i 3, 5 which aUow uniquely to construct the Ai{s^ u). 



Bi{s, u) =ei/^^|2(s - m^)Gi{s, t, r) - 4^4(5, t, r) - 2rG5(s, t, r) 

2 

- 2tGQ{s, t, e) + 2^4(5, u, e) + [274(5, k") - 276(s, e) - ri(s, k") 



+ r2(s, + 2r4(s, k"") - 2Vq{s, e)] + ^ [-71 (s, /c^) + 72(s, A;^) + 275(5, A;^) 

- 276(s, k") - ri(s, + r2(s, k") + 2r5(s, /c") - 2r6(s, k")] I + (s ^ , 



52 (s, u) =e^3^|2fi(t) + (s - m^) [2^2(5, t, A;^) _ 2^3(5, t, A;^) + f]2(s, A;^) 

- 2^6(5, u, P)] + 4^4(5, t, A;2) + 2^2^5(5, t, A;2) + 2(2m2 - s - u)G(^{s, t, k'^) 

+ 2(M2 - t - k'^)G7{s, t, k^) + 2tG9(s, t, A;^) + 2{u-s- 4m^ + - k'^)G8{s, t, k^) 

+ 2(s - - [2^10(5, t, A;2) + l]io(s, u, A;^)] + 2(M2 - Sm^ - s)l]8(s, A;^) 

^2 

- 2f]4 (s, A;2) + 272 (s,A:2) + ^[71(5, A;^) -72(5, A;2) -275(5, A;2) 

+ 276(s, A;2) + ri(s, - r2(s, A;^) - 2r5(s, k^) + 2r6(s, A;^)] | + (s ^ , 



5i(s, «)+2S4(s, u) = ey^^|(s - m^) [2^1(5, t, k^) + 2^2(5, t, A;^) - 2^3(5, t, k^) 



- 405(5, t, A;^) + ni(s, -u, A;^) - 03(5, u, k^) - 205(5, w, k^)\ + 2(M^ - 3m^ - s) 

■ [2Gq{s, t, k^) - 2^10(5, t, A;2) + 06(s, A;^) - Oio(s, ti, A;^)] + 2(2s - 2m^ - Ml) 

■ [2^7(5, t, A;2) + 07(5, w, A;^)] + 2(771^ + - s) [2^9(5, t, A;^) + ^9(5, A;^)] 

+ - '-^^'^'^ + ^2^ + + 2^^^"' 

- 2r6(s, A;2) + 72(5, A;2) + 275(5, A;^) - 275(5, A;^)] | - (s ^ , 

^6(5, U) = —gl^{T^{s, e) - T2{S, k") + 2r6(5, k^) + 270(5, k^) 

tirt 

- 8m'^Gs{s, t, k"^) - 4m^08(5, w, A;^)} + (5 ^ tt) , 
29 



Bs{s, u) =-glN{ri{s, k") - T2{s, k") + 2T^{s, k") + 276(s, k") 
m 

— 8m?Gio{s, t, A;^) — 4m^f2io(s, u, A;^)} — (s <-> tt) . 



(^.7) 



A.2. THE (-) AMPLITUDES 

Electric and magnetic Born terms: 

A[-\s,u) = -Bi-\s,u) = ^g^^f^-^ L^) F^{k^) , Ai-^{s,u) = 0, 

Bi-\s,u) + 2Bi-\s,u) = J^^F;{e) + '-^ [FUe) - F^ie)] , (-4.8) 

Here, F^{k'^) = Ff{k'^) — F^{k'^) denotes the nucleon isovector Dirac form factor and 

is the one loop representation of the pion charge form factor. It is important to note 
that the Born amplitude is more than just the tree amplitude (with form factors equal to 
unity) multiplied by the appropriate form factors generated by the loops. It is known that 
such a simple multiplication prescription violates gauge invariance, unless the (unphysical) 
constraint F^{k'^) = F^{k^) is imposed (or one adds some additional seagull terms [2,16]). 
Within the framework of CHPT the form factors generated by loops and counterterms are 
automatically included in a gauge invariant fashion. 

472!5,6(«:«) = 0, 

4 ^(^'^) = -^^^-(7^ + ^73;;^)' (AlO) 

A'f\s,u) = - ^^SKy[ — ^ —^1 ■ 

4m \ s — 777.^ u — J 

where 5kv = SKp — dKn is adjusted to reproduce the empirical value of the nucleon isovector 
anomalous magnetic moment, Ky = 3.71. 
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For class I and II the following representation of the (-) amplitudes is most economic, 

a[^\s,u) = ai{s, k^) - ai{u,k'^) - B^~\s,u), 



(-) k"^ P 

s — u — 

B^~\s, u) + 2B^^\s, u) = 6i4(s, k^) + hi^iu, k^), A^f\s, u) = 03(5, fc^) + asiu, k^), 

A^^\s,u) = a3(s, k^) - as{u, k^), A'q~\s,u) = ae{s, k^) + aQ{u, k'^) . 

(All) 

Class I: 

ai(s,/e^) = ^^(s-m^)[76(s,/c^) -74(s,/c^)], 

TT 

614(5, k') ^^{{s + 3m2)[276(., k') - 72(^, k')] + {s - m')Ms, k') - 2^,{s, k')]}, 

TV 

03(5, fc^) = ^^|^m74(s,fc^), 

TT 

aeis, k^) = ^y^m[2^6{s, k^) - 72(5, k^)] . 

{AA2) 

Class II: 

3 

ai(s,A;2) = ^{(s-m2)[874(s,A;2)_ 876(s,A;2)-ri(s,A;2) + r2(s,A;^) 

+ 2r4(s, k^) - 2Vq{s, k^)] + 8m2[474(s, k^) + r4(s, A;^)]}, 

614(5, '^') = ^{(^ + 3m')[875(«, k^) - 471 (^, k^) - ^i{s. k') + T^is, k^) + 2V^{s, A:^)] 

+ ^' + lQ^^' + 5^' [4^^(,, fc2) _ 8^^(^^ ^2) _ 2r,(,, e)] + 4m2[-875(fc^) 
s — 

+ 471 (A;2) + fi(P) _ 2t^{k^) + 8m2(-27^(A;2) + 47^(A;2) + r(i(P))]}, 

3 

03(5, /c") = ^{[876(s, fc") - 1674(5, k") - 4r4(s, /c") + ri(s, /c^) 

-r2(s,A;2) + 2r6(s,/c2)-^^[r4(s,A;2) + 474(5, 

3 

a6(s, /c') = ^{875(5, k") - 471 (s, /c2) _ ri(s, /c2) + r2(s, A;2) + 2T^{s, k") 

+ ^(^ + ^') [27,(., e) - 476(., _ r6(., e)] - 8Uk^) + 47i(/c2) 
s — 

+ fi(A;2) _ 2f 5(/c2) + 8m2[-27^(A;2) + 47^(A;2) + r(j(A;2)]} . 

(A.13) 
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Here, the prime on the loop functions denotes the partial derivative with respect to s 
evaluated at s = m^. Note, that we have 6i4(m^, k"^) = aQ^m'^, k'^) = such that i?2 (■§, u) 
has no more pole at s = or u = m^. The residue at these poles is proportional to 
the isovector Dirac form factor, already contained in the electric Born term. Similarly, 
4"^(m2,m2) = 0. 

Class III: 



Zi \ S TTt 



+ rG5(s, t, r) + (s - u)Gq{s, t, k^) + (M^ - r - t)G7(s, t, k^) + tG9(s, t, r) 

w? 4 
4m^ 



2 7 2 

+ (W - t)G8(s, t, k") + 2^ k-") - V2{s, k') + 2Teis, e)] 



Z \ S TTt 

- G3(s, t, /c^)] - 2G4(s, t, /c^) - /c2G5(s, t, k^) + {s + u- 2m'^)G6{s, t, A;^) 
+ (t + - M^)G7{s, t, k^) + {s-u + 4m^ - M ^ + k^)G8{s, t, k^) 

k^ 

- tGg{s, t, + 2{m^ + Ml- s)G^ois, t, k'') + [-Ti{s, k") 

+ r2(s, k") + 2r5(s, k") - 2r6(s, k")] | - (s ^ , 

+ 2G5(s, t, /c^)] + 2{s -Ml + 3m2) [Gels, t, A;^) - Gio(s, t, /c^)] 
+ 2(2m2 + M| - 2s)G7(s, t, A;^) + 2(s - - MI)G<){s, t, k^) 

2 

e — Tfl 

+ [ri(s, e) - T2{s, k") - 2T^{s, e) + 2Tq{s, k")] 

- W + [(^ - 4m2)f + (4m2 - M^)? ^(M^)] | + (s ^ , 



tx) = —gl^[-T^{s, e) + r2(s, A;2) - 2Tq{s, k^) + Sm^Gsis, t, A;^)} - {s ^ u) 
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Am 



B^s = ^glNi-'^iis,^^) + ^2{s,k^) - 2re{s,k^) + 8m^G^o{s,t,k^)} + {s ^ u) . 

{A.U) 

It is important to note that in the chiral hmit = one has the relation 



2mt 



l-—GA{t) 
9A 



(yl.l5) 



with GA{t) the one loop expression for the nucleon axial form factor. This relation contains 
the basic current algebra statement, that the nucleon axial form factor enters the pion 
electroproduction amplitudes. Since in general current algebra statements are exact only 
in the chiral limit we can not expect eq.(A.15) to hold also for finite pion mass. In that 
case corrections beyond current algebra show up which are not constrained a priori. 



APPENDIX B: LOOP FUNCTIONS AND THEIR IMAGINARY PARTS 

Here, we will heavily borrow from appendix B of ref.[l] where all the loop function 
occuring in photoproduction have been defined and given in terms of their Feynman pa- 
rameter representations. We will write down here only the relevant extensions to A;^ 7^ 0. 
The one loop expressions of the form factors involve the loop function 



J 22 (^^) — 



1 ffc^ 



— + / dx[M^ + k^x{x - 1)] In 



327r2\ 6 jQ 

Furthermore, we have in our electroproduction amplitudes 

1 



(S.l) 



li{s,k ) = 



, / dx f'dy^-^^y^^ 

167r^ Jo Jo hy{x,y;s,k^) 



{i ^ 3), 



{B.2) 



with hj{x,y;s,k'^) = M^(l - y) + m?y'^ + {s - m?)xy{y - 1) + /c^(l -yYx{x- 1) and the 
Pi{x,y) are the same as in ref.[l]. An other set of new functions is: 



ri{s,k^)^ 



167r2 



dx / dy— 

Jo hr[x,y;s,k^) 
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{^ ^ 3), 



(S.3) 



with hr{x, y; s, /c^) — M^{1 — y) +w?y'^ + (■^ — w?)xy{y — 1) + k'^y'^x{x — 1) and the qi{x, y) 
are the same as in ref.[l]. For the nucleon form factors their values at s = occur which 
we denote by a tilde, 

= li{m\ e) , fiik") = Ti{m^, k") , {i ^ 3), 
73 ) = ^^T^ / dy{l-y)lii-^ —- , , . 



hr{x,y;m'^,0) 



and an overbar means subtraction at fc^ = 0. 

The box graphs give rise to the G-functions, 

,1.2^ 1 f\ [\ f' , {l-x)yHi{x,y,z) 

Gi{s,t,k) = —^ dx dy dz—^ TT^ ^ ^ 4), 

167r^ 7o Jo Jo hfj{x,y,z;s,t,k^) 

1 (x-l)y^ 

G4(^,t,/c')= 32"-^ / dx dy dz y . t k'^V 
<J^~ Jq Jq Jq nGyX,y, z, s,i, k ) 

with 

hG{x,y,z; s,t, /c^) =M^[1 - y + x{x - l)y'^z] + m'^y'^ + {s - m^)xy{y - 1) 
+ t{l - xfy^z{z - 1) + k'^x{x - l)y2(l - z) 

and the fi-functions, 

.2^_ 1 [\.. f\i.M-y)yn{x,y,z) 



n^is,u,k) = ^,l^ dxj^ dyj^ (^^4), 

1 r... iy-^)y 



J]4(s,w, A;^) = J dx j dy J dz 



Jo hn{x,y,z;s,u,k'^) 



with 



(S.5) 



(S.6) 



(S.7) 



(5.8) 



hn{x, y, z; s, u, /c^) =M|[1 - y + x(x - l)y^] + m^y^ + - m^){l - x)y{y - l)z 

+ (u - m'^)xy{l -y){z-l) + k'^{l - yfz{z - 1) . 

The polynomials bi{x,y, z) and ri{x,y,z) are the same as in ref.[l]. 

Finally, we give the expressions for the imaginary parts of the new A;^-dependent scalar 
loop functions (those with i = 0). The imaginary part of the vector and tensor functions are 
obtained in the standard way. One has to write down appropriate linear relations among 
these functions. The solution of a subset of these linear equations with maximal rank 
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gives rise to explicit formulae for the imaginary parts of the vector and tensor functions 
expressed through those of the scalar ones. 

The quadratic polynomial \{x,y,z) — x'^ + y'^ + z'^ — 2xy — 2xz — 2yz is called the 
Kallen or triangle function. Let us abbreviate = X{s,m?, M^) and = X{s,m?,k''). 
Using the Cutkosky cutting rules we obtain the following imaginary parts: 



T r r = ^ 1 (g + ^ - M^){s -m' + k') - 2sk' + VA^Afc 
167rVX^''{s + m^-M^){s-m^+k^)-2sk^-VX;X;' 



ImGn s, t,k ) = , . m , , 

327rV5^ - AgCg Ag^Bg-^BI-AgCg 

Ag{s, t, k') = m^s'^ - s(2m^ + k^M^) + m^ + k'^M^ik'^ + - 3rr?) , 
Bg{s, t, e) = -^A, + \mI{MI - e){Ml + 2e - Sm^ - s) , 
CG{s,t,e)=tX^ + M^{M^-kY, 



(S.9) 



Imnois, u, k') = m ^ii±^^i±v^^S , 

327r/B^^^A^ A^+Ba-sjBl-A^Cn 

An{s, u, /c^) = M^(s - m'f + k'M^{M^ - 3m^ - s) + kV , 

Bn{s,u,k') = -|a. + ^(m^ - M^) + s{-m' + 2m'M^ - ^M^ - m'k' + ^M^k') 



1 .6 7 



+ -m" - -m^M^ + 2m' + k'{m^ + -m'M^ - -M^) - k'^m' , 
2 2 2 2 

Cn(s,w,/c2) = wA^ + 2s(M2 - m'){Ml - k') + {2m' - - k'){3m'M^ - 2M^ - k'm') 



The formulae given above hold for fc^ < 0, t < 0, u < {m — Mj^)' and s > (m + 
M^)2. The physical region is completely contained within this domain of the Mandelstam 
plane. 

APPENDIX C: CALCULATION OF CURRENT MATRIX ELEMENTS 
FROM MULTIPOLE AMPLITUDES 

To calculate polarization observables using the formulas of Ref.9, we need to construct 
current matrix elements in the center-of-mass frame of the final ttN system. The photon 
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momentum defines the z-axis and hence we have /c^ = {ko,0,0, k). The amphtude of the 
7*A'" —>■ ttN process can then be written as 



< rus' I e^{X)J^ I >= < rus' \ Oi \ rus > fi{k^, W, x) {C.l) 

i=l,6 



with 





— ia 


■ ex 




02 


= a ■ 


qa- 


(k X ex) 


03 


= ia 


■ kq 


• e\ 


04 


= ia 


■ QQ 




05 


= ia 


■kk 




Oe 


= ia 


■kk 





{C.2) 



Here we have defined q as the pion momentum in ttN center-of-mass frame, q — q/\q\ and 
k = k/\k\. The spherical unit vectors are defined as e±i = ^l/\/2{x±iy) and eo = The 
photon polarization vectors for A = ±1, in Eq.(C.l) are defined as 



e'^(±l) = (0,e-'±i) 



1 



:{\k\,koeo) 



(as) 



where k"^ = k^ — k"^ < 0. The coefficients /i(/c^, W,x) are only functions of ttN invariant 
mass W , four-momentum transfer square A;^ and x = k ■ q. In terms of the multipole 
amplitudes, we have (all multipoles are functions of k'^ and W) 



h{k^W,x 
h{k\W,x 
h{k\W,x 
f^{k\W,x 
h{k\W,x 
h{k\W,x 



= J2iEi+Pl+,{x) + Ei_PU{x) + Mi+Pl+,{x) + l)Mi_PU{x)) 

I 

I 

= Y,{Ei+Pl'^^{x) + Ei.PlUix) - Mi+Pl'^^{x) - Mi.P['_Ax)) 
I 

= ^{-Ei+P;'{x) - Ei_Pl'{x) + Mi^Pl'ix) - Mi_P;'{x)) 
I 

= j2i-ii+m+pii^)+isi-Piix)) 

I 

= + m+pi+ii^) - isi-puix)) 



{CA) 
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By using Eq.(C.3), it is easy to see that 

< nis' \ Jx \ rris > = -1/V2(< nis' \ e^{+l)J^ \ tHs > - < nis' \ e^(-l)J^ | nis >) ,^ , 

(CSa) 

< nig' \ Jy \ nis > — i/V2{< nis' \ | > + < nis' \ ^^{—^)Jn \ 'rris > 

By using the current conservation condition koJo = k ■ J = kJg, Eq.(C.l) for A = leads 
to 

<ms' \ J^\ms> = < ms' \ e^(0)J^ \ nig > 

Kq 

The matrix elements defined in Eq.(C.5) are needed to calculated various polarization 
observables defined in Ref.9. 

APPENDIX D: MODIFICATION OF INVARIANT AMPLITUDES IN THE 
CASE OF UNEQUAL CHARGED AND NEUTRAL PION MASSES. 

In ref.[7] is was observed that isospin breaking effects due to the difference of the 
neutral and charged pion masses turn out to be quite sizeable in the case of neutral pion 
photopro duct ion close to threshold. From this experience one expects that this isospin 
breaking effect will also be relevant for 7r° electroproduction from protons. Therefore we 
will display here the relevant modifications of the invariant amplitudes Ai{s,u) for the 
reaction 7*p Tc^p which arise if we do not set equal the neutral and charged pion mass 
as it was done in the chapter IV. Notice that in contrast to the photoproduction case, one 
has to consider diagrams from all three gauge invariant classes. 

Class I and II: 

The expressions given in eqs. (4.5,4.6,4.7) hold still for the invariant amplitudes. There 
is only some change in the loop functions ^i{s,k^) and Ti{s,k'^) of the following form. 

7i(s, /c^) has to be evaluated with M^r = and is denoted 7i+(s, A;^) whereas Ti{s,k'^) 
has to be calculated with = M^o denoted by rjo(s, k^). 

Electric and magnetic Born terms: 

Here, the proton Dirac form factor Ff{k'^) in eq.(3.12) gets modified. The 'yi{k'^) and 

~ ~ — 7T7T 

Ti^k^) are replaced by 7j_|_(/c^) and Fio^k'^), respectively. Furthermore, J22 (^^) of eq.(B.l) 
has to be evaluated with M^r = Mt^+. Of course, now the counter term contribution 
proportional to 6g(A) + 6g(A) has to be readjusted in order to reproduce the empirical 
proton mean square charge radius. In a similar fashion the value of the magnetic moment 
counter term SKp is obtained from 

5Kp = Kp + 8£^^J^m^{^4o(m^0) -2 74+(m^0)} = 1.291 (DA) 

Class III: 

Most of the terms in the Bi{s,u) can be taken over from eq.(4.8) by simply re- 
placing the respective loop functions by 7j_|_(s, A;^), rjo(s, /c^), rio(t), Gio{s,t, k"^) and 
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u, k"^). In all kinematical prefactors like — t — k'^, M^r reads M^o. The modified 
loop functions need some explanation since they involve both the neutral and charged 
pion. The proper denominator for their Feynman parameter representation is now 

hn+ {x, y, z; s, u, k^) =M^+ (1 - y) + M^ox{x - + mY + (s - m^) (1 - x)y{y -\)z 

+ - m^)xy{l - y){z - 1) + k^{l - yfz{z - 1) . 

(D.2) 



The amplitudes i?i^2,3(-s, w) involve the combination 



/Lf2 

rr (^) - ^^r^^(.) - vi^{s) . (i?.3) 

5 777' 



It gets replaced by 



2M^ 2 



where the Feynman parameter representation of the loop functions reads 

^+^^^ IGTr^io ''io ^M2+(l-y) + M2oxy2(x-l) + mV + (^-™')^y(y-l) 
with ^0 = I5 ^1 = 1 — 2/5 5'2 = 1 — a^J/ and rfg^(s) is obtained if is set equal M^o. 
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FIGURE CAPTIONS 



Fig.l. The functions Si^2,3,4(p) defined in eq.(5.2) with p = —k'^/M^. The sohd, dotted, 
dashed and dash-dotted hnes refer to Si, S2, S3 and S4, respectively. 

Fig.2. The functions $i,2(/c^) defined in eq.(5.5). The sohd and dotted hnes refer to $1 and 
$2, in order. 

Fig.3. The functions Si^2(po) defined in eq.(5.9) with po = —k'^/M^o- The sohd and dotted 
hnes refer to Si and S2, respectively. 

Fig. 4. Comparison of the theoretical predictions with the recent NIKHEF data for ^(7*, 7T^)p 
[3]. The solid circles, open circles and open squares refer to the one-loop CHPT, tree 
and PV Born term calculation, in order. 

Fig.5 The longitudinal So+ multipole amplitudes for p{j*, TT~^)n, n(7*, tt~)p and ^(7*, TT^)p. 
The solid, dotted and dashed lines refer to the one-loop CHPT, the tree and the PV 
Born calculation, in order. 

Fig.6 The transversal £'0+ multipole amplitudes for ^(7*, 7r"'')n, 77,(7*, 7r~)p ^nd p(7*,7r°)p. 

Fig. 7 fc^-dependence of ctt and ctl for ^(7*, 7r'^)n, n(7*, 7T~)p and ^(7*, TT^)p. 

Fig. 8 Angular distributions of the transverse cross sections for p{j*,TT~^)n, n{j* ,tt~)p and 
p(7*, 7r°)p. For notations see Fig.5. 

Fig. 9 Same as Fig. 8 except for the longitudinal cross sections. 

Fig. 10 Same as Fig. 8 except for the interference cross sections. 

Fig. 11 Same as Fig. 8 except for the polarization cross sections. 

Fig. 12 Dependence of the axial cut-off mass for ^(7*, 7r"'")?i. We show the angular distri- 
butions dar/dQ, dai/dQ, daj/dQ and daj/dQ at /c^ = -0.04 GeV^ for Ma = 
0.96, 1.06, 1.16 GeV corresponding to the dotted, solid and dashed lines, respectively. 

Fig. 13 The S-wave cross section ao for p{'y*,n^)p including isospin-breaking for W = 1074 
MeV and e = 0.58. The solid, dotted and dashed lines represent the isosymmetric 
CHPT, CHPT with isospin-breaking and PV predictions, in order. 

Fig. 14 Angular distributions for p{j*,TT^)p including isospin-breaking. For notations see 
Fig.l3. 

Fig. 15 /c^-dependence of (Jt and (Jl for p{'y*,TT^)p including isospin-breaking. For notations 
see Fig. 13. 
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TABLES 





-0.001 GeV2 


-0.04 GeV2 


-0.06 GeV2 




( 25.77, 0.2663) 


( 15.00, 0.5830) 


( 9.899, 0.6602) 




(-18.44, 0.1890) 


( -21.32, 0.2296) 


(-22.58, 0.2425) 


Ei^ 


(0.7625, 2.022 10"^) 


(0.6366, -1.169 10-6) 


(0.5707, -7.883 10"^) 


Si+ 


(-0.4098, 2.631 10"^) 


(-0.3335, 7.441 10"^) 


(-0.2849, 9.871 10"^) 


Si- 


(-3.373, 9.221 10-^) 


(-4.831, 0.01405) 


(-5.152, 0.01523) 


Mi+ 


(-1.037, -5.792 10"^) 


(-1.135, -9.365 10"^) 


(-1.168, -1.058 10"^) 


Mi_ 


(0.2250, -6.891 10"^) 


(-0.8534, 2.085 10-^) 


(-1.271, 2.989 10-^) 



Table la: Various multipoles for the reaction — > Tr+n in units of 10 ^/M^+. 





-0.001 GeV^ 


-0.04 GeV^ 


-0.06 GeV^ 




(-29.99, -0.3543) 


(-16.96, -0.6121) 


(-10.86, -0.6657) 




(22.98, -0.08447) 


(26.55, -0.08553) 


(27.24, -0.08384) 


Ei+ 


(-0.7534, -4.067 lO'^^) 


(-0.6263, 3.178 10"^) 


(-0.5611, 5.027 10-9) 


Si+ 


(0.4010, -1.952 10-^0) 


(0.3174, -7.329 10"^°) 


(0.2691, -1.377 10-9) 


Si- 


(3.336, -8.639 10-^) 


(4.778, -0.01237) 


(5.137, -0.01300) 


Mi+ 


(1.242, 3.158 10-") 


(1.350, -3.199 10-9) 


(1.341, -5.050 10-9) 


Mi_ 


(-0.5396, 1.667 lO'^) 


(0.5531, 3.483 10"^) 


(1.059, -1.967 10-4) 



Table lb: Various multipoles for the reaction 7*n —>■ n p in units of 10 ^/M^+. 
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A;2 


-0.001 GeV2 


-0.04 GeV2 


-0.06 GeV2 




(-1.304, 0.5103) 


(1.625, 0.9461) 


(2.668, 1.045) 


c 

'-'0+ 


/o inn r\ i o^F;^ 

(o.iyy, U.ioUo J 


[o. ibo, U.ziUz J 


(o.ooD, U.zi (b) 


Ei+ 


(6.989 10-3, -2.634 10"^) 


(1.609 10-3, -4.267 10"^) 


(-2.648 10-3, -4.385 10"^) 


Si+ 


(-0.01311, -7.260 10-^) 


(-0.03610, -7.342 10"^) 


(-0.04627, -7.053 10"^) 


Si- 


(0.09471, 0.01126) 


(0.2800, 0.01760) 


(0.3872, 0.01912) 


Mi+ 


(0.1162, 2.884 10"^) 


(0.04723, 9.844 10"^) 


(-0.02483, 6.400 10"^) 


Mi_ 


(-1.521, -1.233 10-^) 


(-2.303, 2.877 10-^) 


(-2.520, 4.016 10-3) 



Table Ic: Various multipoles for the reaction 7*7? — > 7r°7? in units of 10-3/M^+. 



MA[GeV] 


1.06 


1.16 


0.96 


Re£;o+[10-3/M^+] 


15.00 


15.98 


13.71 


Re>So+[10-3/M^+] 


-21.32 


-20.42 


-22.51 




10.81 


12.17 


9.04 


cri,[//b] 


16.87 


15.51 


18.70 


(ctt + crL)[lM 


27.68 


27.68 


27.68 


(Tt/ctl 


0.64 


0.79 


0.48 



Table 2a: Dependence on the axial cut-off mass Ma of various observables for 
the reaction 7*^ — > Ti+n. 
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MA[GeV] 


1.06 


1.16 


0.96 


ReEo+[10-yM^+] 


-16.96 


-17.94 


-15.67 


Re,So+[10-7M^+] 


26.55 


25.65 


27.74 




13.70 


15.27 


11.74 




25.58 


23.95 


27.88 


(ctt + (Ti,)[//b] 


39.28 


39.21 


39.62 




0.54 


0.64 


0.42 



Table 2b: Dependence on the axial cut-off mass Ma of various observables for 
the reaction 7*n —>■ TT~p. 
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